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ABSTRACT 

We study in detail the mid-infrared Spitzer-IRS spectrum of M 87 in the range 5 to 20 fim. Thanks to 
the high sensitivity of our Spitzer-IRS spectra we can disentangle the stellar and nuclear components of 
this active galaxy. To this end we have properly subtracted from the M 87 spectrum, the contribution 
of the underlying stellar continuum, derived from passive Virgo galaxies in our sample. The residual 
is a clear power-law, without any additional thermal component, with a zero point consistent with 
that obtained by high spatial resolution, ground based observations. The residual is independent of 
the adopted passive template. This indicates that the lO/im silicate emission shown in spectra of 
M 87 can be entirely accounted for by the underlying old stellar population, leaving little room for 
a possible torus contribution. The MIR power-law has a slope a ^ 0.77-0.82 (Sj^ oc iy~°'), consistent 
with optically thin synchrotron emission. 

Subject headings: - galaxies: individual (M 87) - galaxies: elliptical and lenticular, cD - galaxies: 
active - galaxies: jets - Infrared: galaxies 



1. INTRODUCTION 

iBressan et al.l (|2006D presented high signal-to-noise ra- 
tio Spitzer Infrared Spectrograph (IRS) observations of 
18 Virgo early-type galaxies (ETGs hereafter), selected 
from those that define the color-magnitude relation of 
the cluster, with the aim of detecting the silicate emis- 
sion of their dusty, mass-losing evolved stars. They found 
that 13 ETGs, i.e. 76% of the sample, show a pronounced 
broad silicate feature that is spatially extended and likely 
of stellar origin, i n agreement with m odel predictions of 
old populations (jBressan et al.l Il998l ) . The absence of 
any other feature in their spectra suggested that these 
galaxies could be considered as the prototype of passively 
evolving ETGs. In contrast, the IRS spectrum of the 
other four ETGs presented various levels of activity. 

Among these, NGC 4486 (M 87), the well-known gi- 
ant elliptical located at the cluster centre, stands out. 
Its mid-infrared (MIR hereafter) spectrum, besides nar- 
row atomic emission lines, shows an excess of emission at 
longer wavelengths with respect to passive ETGs. The 
silicate emission at lO/xm is also evident in the spec- 
trum, and since M 87 hosts a strong central radio source 
(Virgo A) and a jet, visible from radio to X-ray wave- 
lengths, it could be an indication of possible torus emis- 
sion, as already found in other AGNs (Siebenmorgen et 
al. 2005, Shi et al. 2006, Perlman et al. 2007). 

The importance of the torus in the context of the AGN 
energetics makes the debate about the origin of the nu- 
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clear MIR emission in M 87 very intriguing. Early, high 
spatia l resolution, ground-based MIR observations o f 
M 87 (jPerlman et al.ll2001l : IWhvsong fc Antonuccill2004f) 



have shown that the emission from the nucleus is consis- 
tent with synchrotron emission. Both investigations have 
excluded the possibility that M 87 hosts a quasar-like 
nucleus similar to that at the center of Cen A. However, 
neither could exclude the po s sibilit y of a low-luminosity 
AGN. Indeed [Perlman et aP ()2001f ) detected a faint, ex- 
tended, circum-nuclear component that could originate 
in a torus contributing '^7% of the nuclear luminosity. 
Alternatively, they suggested that this extended com- 
ponent could be stellar emission fr om red giant stars. 
More recently, iPerlman et al.l (|2007f ) reanalyzed this is- 
sue by means of high resolution ground based imaging 
and Spitzer spectroscopy. They confirmed the consis- 
tency with synchrotron emission for the MIR fluxes of the 
knot A/B complex in the jet. However, they measured 
a spectral index aiji ~0.41, much lower than before and 
different from that found in the nearby structures (e.g. 
IShi et ai]|2007( ). They still noticed deviations from the 
pure synchrotron emission and speculated that, at least 
in the shorter wavelength range, this could be due either 
to contribution from M and K stars or contamination 
from the low frequency tail of the emission from knot 
HST-1, not resolved by Spitzer but also not apparent in 
the 11.7^m SUBARU image of the nucleus. 

We revisit here the nature of the MIR emission of M 87 
exploiting our high signal-to-noise Spitzer spectra of pas- 
sively evolving ETGs. Taking advantage of these data 
we have constructed very accurate MIR templates of the 
stellar populations in these galaxies (Valdes et al. 2009, 
in prep.). 

These templates will be used to disentangle the nuclear 
non-stellar emission of M 87 by a proper subtraction pro- 
cess. 

In section [2] we summarize our Spitzer IRS observa- 
tions and data extraction as well as all supplementary 
data needed to perform the subtraction. In section [3] we 
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Fig. 1. — IRS SL (nearly horizontal) and LL (nearly vertical) 
slits superimposed on the MIPS image of M 87 at A=24/im. The 
two additional panels on the right side show a portion of the LL2 
(bottom) and SLl (top) 2-D spectral images. Notice that the for- 
mer does clearly show both the spectra of the nucleus and that 
of the nearest bright knots, whereas the SL slit is not affected by 
these components. The extraction width of the LL spectra does 
not include the contribution from these knots. 



describe the subtraction of the spectral energy distribu- 
tion (SED) of a passive ETG from the IRS spectrum of 
M 87 while, in section [4] we discuss the properties of the 
residual component. Our conclusions are presented in 
section [51 

2. OBSERVATIONS 

In this section we describe how we have obtained the 
NIR/MIR SEDs of M 87, and of three other elliptical 
galaxies, NGC 4365, NGC 4382 and NGC 4636. These 
were selected to be "passive", insofar as they show no 
evidence of either recent past star formation or AGN ac- 
tivity. They will be used as templates for the stellar com- 
ponent in our attempt to isolate the non stellar emission 
of the nucleus of M 87. 

2.1. Spitzer IRS 

Standard staring mode short (SLl and SL2) and long 
(LL2) low resolution IRS spectral observations of M 87 
and other selected Virgo cluster ETGs were obtained 
during the first Spitzer General Observer Cycle on 2005 
June 1 as part of program ID 3419 (PI A. Bressan). The 
spectra were extracted in a fixed aperture (3.6" x 18" 
for SL and 10.2" x 10.4" for LL), using a purpose built 
program that allows the reconstru ction of the intrinsic 
galaxy profile (jBressan et al.|[2006l ). 

Position angles of the SL (3.6" wide) and LL (10.2" 
wide) IRS slits in the case of M 87, are shown superim- 
posed to the MIPS 24/im image, in Figure [T] At these 
wavelengths the emission of NGC 4486 appears quite 
complex even at the low resolution of MIPS 24/im. Perl- 
man et al. (2001) have provided a thorough description 
of the MIR spatial emission of the central regions of M87. 
Their 10.8 //m image shows the "nucleus" and the D, F 
and A, B and C knots, whose fluxes are reported in their 
Table 1. To check which of these external components 
can affect our spectra we show the 2-D spectral images of 
of the LL2 and SLl segments in the bottom and top pan- 
els of Figure [U respectively. In the LL2 slit we can easily 
recognize both the spectra of the nucleus and that of the 
nearest bright A, B and C knots, whereas the SLl slit 



is clearly not affected by these three components. Thus, 
the only extra-nuclear component that may contribute 
to our SL spectrum is knot D. However this feature is 
located 3" from the nucleus, has low level flux (7% of the 
nucleus, Perlman et al. 2001) and thus is likely to pro- 
vide only a negligible contribution to the SL spectrum. 
The other components, located at greater distances, are 
well outside the SL slit. On the other hand our LL slit 
contains the strong knots A, B and C that are clearly vis- 
ible in the spectral images. In this case the extraction is 
limited to a width of 10.4" in order to avoid these strong 
knots that are not present in SL. As we will see later our 
results are independent of the inclusion or exclusion of 
the LL segment from our analysis. 

The estimate of the flux emitted by an extended source 
within a fixed slit aperture involves the deconvolution of 
the received flux with the PSF of the instrument. This 
correction is important to obtain the shape of the in- 
trinsic spectral energy distribution (SED) of the galaxy, 
because from the slit loss correction we have estimated 
that for a point source the losses due to the PSF amount 
to ~20% at 5fim and to ~40% at 15/im. Conversely, 
a uniform source will suffer no net losses. In order to 
recover the intrinsic SED we have convolved a surface 
brightness profile model with the instrumental PSF, and 
we have simulated the corresponding linear profile along 
the slits, taking into account the relative position angles 
of the slits and the galaxy. The adopted intrinsic model 
profile is a wavelength dependent two dimensional mod- 
ified King law (Elson at al. 1987): 

-7/2 
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(1) 



(i?c X b/a)\ 

where X and Y are the coordinates along the major and 
minor axis of the galaxies, b/a is the axial ratio taken 
from the literature. Iq, i?c and 7 are free parameters 
that are functions of the wavelength and are obtained 
by fitting the observations with the simulated profile. In 
order to get an accurate determination of the param- 
eters of the profiles several wavelength bins have been 
co-added. This procedure allows us (i) to recog- 
nize v^rhether a particular feature is resolved or 
not, (ii) to reconstruct the intrinsic profile of the 
source and (iii) to account for the vifavelength de- 
pendent PSF losses, in a more general way than 
for point like sources. 

Examples of the fitting procedure for NGC 4486 at two 
selected wavelengths are shown in Figl2l The dashed 
lines show the reconstructed one dimensional intrinsic 
spatial emission profiles (in electrons/s) sampled along 
the slit, while the solid lines represent the corresponding 
one dimensional PSF convolved models that are fitted 
to the data. The error bars in the figure show l-a un- 
certainties that were computed taking into account the 
Poissonian noise from the source and the background, 
the readout noise and dark current noise, as described 
in the Spitzer Observer Manual. The spatial emission 
parameters obtained for all the galaxies used in this pa- 
per are shown in Table 1. At short IRS wavelengths the 
spatial emission of NGC 4486 appears complex. It is a 
superposition of an unresolved central component and 
an extended component. Indeed, in fitting only a sin- 
gle component to the data we find an unresolved central 
source with _Rc=0.r'. Nonetheless, there is an evident 
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TABLE 1 

MIR Spatial emission parameters of selected early type 

GALAXIES 











14.0;im 


GALAXY 


Rc 


7 


Rc 


7 


NGC 4486 


0.1 


0.5 


0.1 


0.7 


NGC 4365 


3.5 


2.0 


4.0 


2.2 


NGC 4382 


1.3 


1.4 


1.6 


1.3 


NGC 4636 


3.7 


0.8 


3.5 


0.9 



excess over this unresolved component beyond ~ 4", that 
does affect our rectangular extraction, 18" long. The 
excess remains visible up to ^14/im beyond which the 
profile is clearly that of the PSF (i?c=0.1"), and the ex- 
cess in the wings is much less important. From this spa- 
tial analysis we conclude that NGC 4486 is progressively 
more dominated by un unresolved nuclear component as 
the wavelength increases. However, the presence of an 
underlying extended emission, evident at 5 fim, can be 
detected even at 9/zm. Since we have shown that the 
external jet does not contribute to the SL segment, we 
argue that this emission is of stellar origin. 

2.2. Other Observations 

IRAC fluxes were measured from archival "post ba- 
sic calibrated data" images using a customized program 
written in IDL. Before flux measurements, the back- 
ground in each image was determined as the modal sky 
value of the off-source portion of each image. Fluxes 
were then measured in circular apertures centred on each 
galaxy. The error on each flux measurement was cal- 
culated as the quadrature sum of the calibration error 
(taken to be 10%) and the standard deviation of the 
background value. The calibration error is the dominant 
term. 

In order to place M 87 IRAC and IRS data on the 
same flux scale, we measured IRAC 5.8 /im fluxes in the 
circular aperture containing the same flux sampled by 
IRS at the same wavelengths This aperture has a radius 
of 4" and has been selected as the reference aperture 



for any band. Consequently, all broad band fluxes in all 
other galaxies were measured within the same aperture 
and their IRS spectra were normalized to their flux in 
IRAC channel 3. 

In the case of M 87, we have also measured the emission 
at 24/im from the MIPS image, which may serve as a use- 
ful check for our IRS extraction procedure. Since the IRS 
extraction procedure indicates that toward these wave- 
lengths M87 is dominated by the nuclear component, we 
have taken into account the effects of the PSF to ob- 
tain the MIPS 24/xm flux. To this end we have divided 
the flux measured within the nuclear 4" radius aperture 
with the fractional energy of the MIPS 24/im PSF image, 
encircled by the same radius. 

NIR images of the sample galaxies in the J, H and Kg 
bands were downloaded from the All-Sky Data Release 
from the Two Micron All Sky Survey. ^ 

Using these data, in addition to our IRS spectra, we 
have obtained the NIR/MIR SED of M87 and of the 
other three ETGs that will be used as templates for the 
subtraction of the stellar component. 

3. STELLAR COMPONENT SUBTRACTION 

In order to disentangle the nuclear and the extended 
component we resort to the comparison of the high 
signal-to-noise spectral energy distribution (SED) of 
M 87 with that of the passive ellipticals we have also 
observed with Spitzer. Here we describe the process of 
subtraction of the stellar component. 

Before subtraction, all the template galaxy SEDs were 
normaHzed to the flux of M 87 in the 2MASS Ks band. 
This band was selected because it is the nearest band to 
the IRS spectrum with a negligible non-stellar compo- 
nent. The NIR/MIR SED of M 87 and of the template 
ETGs, after this normalization, are shown in Figure [31 

We notice that all SEDs are very similar in the NIR 

NASA/IPAC Infrared Science Archive for NASA Infrared and 
Submillimeter Data chartered to serve calibrated science products 
from NASA's infrared and submillimeter missions. 2MASS Image 
Services provide the users calibrated images that contain photo- 
metric zero point information in their headers and are suitable for 
quantitative photometric measurements. 
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Fig. 3.— Near-IR/MIR SED of M 87 (NGC 4486) and of the 
other three passive galaxies in a circular aperture of 4" radius. The 
passive ETGs have been normalized to the Ks broad band flux of 
M 87. 



TABLE 2 

M 87 CORE MIR POWER-LAW (Si/ (X ) 



Galaxy 


Q 






S3. 6 


S11.7 


Template 






mJy 


mJy 




NGU 4365 


0.766±0.004 


0.756±0.007 


0.06 


8.9 


20.0 


NGC 4382 


0.807±0.005 


0.821±0.007 


0.05 


7.7 


21.3 


NGC 4636 


0.817±0.005 


0.774±0.008 


0.06 


8.1 


22.0 



bands. This similarity, that extends to optical wave- 
lengths, indicates a uniformity of the underlying stellar 
populations. The SEDs of passive ETGs are very similar 
also in the MIR, especially where the integrated spec- 
trum is dominated by the stellar photospheres, i.e. be- 
low 8/xm. At longer wavelengths some differences appear, 
likely caused by small differences in the average proper- 
ties of circu mstellar envelopes of mass loosing AGB stars 
(jBressan e t al. 2006). This uniformity is strength- 
ened by the recent comparison of the surface 
brightness distribution in the Ks and IRS Peak- 
up 16/im bands , of bright ETGs in Virgo and 
Coma clusters ( Clemens et al.l l2009h . NGC 4636 
shows also narrow emission lines but its stellar contin- 
uum is that of a passive elliptical. 

In contrast, M 87 shows a clear MIR excess which 
reaches about one order of magnitude at Ri20/im (see 
also Perlman et al. 2007). An estimate of the quality 
of our IRS extraction process is provided by the inde- 
pendent measure of the 24/im MIPS flux, plotted in the 
same figure. 

In summary, inspection of Figure [3] shows that M87 
stands out in the MIR spectral window. In this region, 
we hypothesize that the stellar populations of M 87 and 
the other passive ETGs are similar, and we subtract the 
spectra of the passive ETGs from that of M 87. 

4. THE RESIDUAL COMPONENT 

To highlight the nature of the MIR excess we sim- 
ply subtract the SEDs of the passive ETGs from that 
of M 87. Figure 2] shows the results obtained using 
the SED of NGC 4382, but similar results are obtained 
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Fig. 4. — MIR residual after subtracting the Ks band normalized 
template SED of NGC 4382 from that of M 87 (NGC 4486), both 
indicated in the figure. The exponent of the power-law fit (dashed 
line, Siy tx is shown. The dot and the triangle mark 

the 11.6 7/^m and 10.8/jm grou nd based nuclear measure- 
ments bv lPerlman et al.l (|2007i^ and jPerlman et al](|2001h . 
respectively. The open square is the M 87 MIPS 24 fim flux. Af- 
ter correction for the NGC 4382 extrapolated emission, the M 87 
MIPS 24 fira flux lowers to the position indicated by the filled 
square. Crosses show another IRS spectrum of M 87 that was 
extracted from the a rchive data (ID 82, PI G.H. Rieke) used by 
IPerlman etal] i200% . 



adopting the other galaxies as templates for the stellar 
emission. The residual, without emission lines that have 
been masked out, is well-fitted with a power-law and does 
not show any evidence of other particular features remi- 
niscent of additional emission components. This is even 
more surprising considering the prominent broad features 
visible in the stellar template and given the relative con- 
tribution of the stellar and non-stellar components. Ta- 
ble 2 summarizes the power-law exponents (S^ oc v^"), 
the fractional contribution of the power-law in the Ks 
band, and its absolute flux at 3.6 and 11.7/im, for the 
three templates adopted. Errors in the exponents are 
only formal errors of the fits and do not include errors 
from the normalizations involved in our procedure. A 
realistic uncertainty is provided by the range of slopes 
obtained with different galaxy templates. 

Our residual flux (Table 2) is consistent with 
th e high resolutio n, ground based nuclear values 
by IPerlman et al.l ((^07,) (S^= 20.8±3.5 mJy at 
11.67Mm, FWHM=0" 55, Figure H red dot) and 
IPerlman e t al.' (2001') (S^ = 16.7 ± 0.9 mJy at 
10.8/im, FWHM= 0'.'46, Figure H red triangle). 
In the same Figure, we have also plotted the observed 
MIPS 24 ^m flux, before (open square) and after (filled 
square) a 15% correction due to the contribution of stars, 
estimated from the extrapolation of the passive elliptical 
template. The MIPS 24 /im flux is also a time indepen- 
dent measure and falls on the same residual power-law 

The power-law has a slope a between 0.77-0.82, de- 
pending on the adopted stellar population template. 
This value is signifi cantly larger than that found by 
IPerlman et all (pOOl (~0.41). They fitted different IRS 
observations (PID82, Rieke), within a lower extraction 
width and using the spectral range 7.5 to 15/im, specif- 
ically selected to avoid the shorter wavelength region 
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where the contamination from the stars could be signif- 
icant. This large discrepancy could be real, for example 
due to a variability of the nuclear component, or could 
simply be due to the different method adopted for the 
analysis. To understand its origin we have r e-analysed 
the IRS data set used bv lPerlman et al.l ()2007l ). with our 
extraction procedure. The resulting spectrum is plot- 
ted in Figure [4] (crosses). We noti ce that our flux value s 
are larger than those published bv lPerlman et al] (|2007f ). 
but this is due to the 2.5 times larger extraction width 
that we have adopted. Since, apart from this, the two 
IRS datasets are very similar, we argue that the differ- 
ent values of the slope arise from the different approach 
adopted. A crude fit to our M 87 SED, in the 7.5 to 
15/im region, gives a slope a ~0.17. But in our case the 
strength of the underlying stellar population is not neg- 
ligible and changes from 70% at 5 /xm to 25% at 15 /im. 
After the subtraction, the slope of the r esidual is a ~ 
0.8. W ithin the smaller aperture used bv lPerlman et al.l 
(l2007t l the contribution of the stellar population is cer- 
tainly less severe than in our case, but not negligible, in 
particular around 7.5 /itm. The lack of this correction is 
most likely the cause of the different slope obtained. 

Finally, we notice that th e slope we der i ved is more 
consistent with the previous iPerlman et"al] ()200lD anal- 
ysis. Our slope is also similar to the value found for jet 
knots A/B and C/G by Shi et al. 2007, and to the slope 
found at radio wavelengths when the entire, extended 
radio source is measured. 

5. CONCLUSIONS 

Exploiting our high signal-to-noise Spitzer IRS spec- 
tra and broad band data from the literature, we have 
built template SEDs of the passive stellar populations 
in ETGs. After a passive stellar component is properly 
subtracted from our M 87 IRS spectrum, the residual, 
almost independent of the adopted template, is a power- 
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law, and lacks any signature of the 10/xm feature. This 
suggests that the 10/xm bump observed in the IRS spec- 
trum of M 87 can be entirely explained in terms of the 
underlying stellar population, typical of all other pas- 
sively evolving Virgo ellipticals (Bressan et al. 2006). 
The residual intensity matches fairly well the high spa- 
tial resolution ground based SUB ARU imc\ea.i measure - 
mer its at lO.SAtni and 11. 67Mni, bv lPerlman et al.l ()200lD 
and lPerlman et all (|2007f ). respectively. Since the contri- 
bution of stars at these wavelengths is not negligible, this 
indicates that our procedure, even if applied to a large 
aperture (r~4" corresponding to «300 pc at 16 Mpc) 
is accurate enough to reveal the emission from the <1" 
core. 

The slope of the power law is strikingly similar to that 
found for the nearby jet knots and for the large aperture 
radio data. 

The above findings confirm, in a new independent way, 
that the inner (<!") core of M 87 is dominated by syn- 
chrotron emission, leaving little room for a torus c ompo- 
nent (iWhvsong fc Antonucc i 2004; Perlm an et al.l l2007'). 

Since the signatures of MIR silicate emission have 
been detected in many AGN observed with Spitzer 
(|Siebenmorgen et al.|[200l IShi et a"l]|2006D . we suggest 
that, in some circumstances, the underlying stellar pop- 
ulation may be an important component to be taken into 
consideration. 
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